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VISCOUS/ POTENTIAL FLOW ABOUT MULTI -ELEMENT r.'.T»-DIMENSIONAL AND 
INFINITE-SPAN SWEPT WINGS: THEORY AND EXPERIMENT 


L. E. Olson* 

Ames Research Center, NASA, and U.S. Army Air Mobility RfiD Laboratory 
Moffett Field, California 
and . 

F. A. Dvorak* 

Analytic Methods, Bellevue, Washington 


Abstract 

The viscous subsonic flow past two-dimensional 
and infinite-span swept multi-component airfoils is 
studied theoretically and experimentally. The com- 
puterized analysis is based on iteratively coupled 
boundary-layer and potential-flow analysis. The 
method, which is restricted to flows with only 
slight separation, gives surface pressure distribu- 
tion, chordwise and spanwise boundary- layer char- 
acteristics, lift, drag, and pitching moment for 
airfoil configurations with up to four elements. 
Merging confluent boundary layers are treated. 
Iheoretical predictions are compared with an exact 
theoretical potential flow solution and with experi- 
mental measures made in the Ames 40- by 80-Foot Wind 
Tunnel for both two-dimensional and infinite-span 
swept wing configurations. Section lift character- 
istics are accurately predicted for zero and moder- 
ate sweep angles where flow separation effects are 
negligible. 


a angle of attack 
6 * sweep angle 

Y * vorticity 

5 = panel angle relative to reference coordinate 

system 

6* * chordwise displacement thickness, ro 

e « momentum thickness, m 

V ■ kinematic viscosity, m^s** 

Subscripts 

e = edge of viscous layer 

exp > experimental value 

f * flap 

i,j = panel indices 

k » component indice 

te ■« trailing edge 

vlt = vortex- lattice theory 

w « wing 

= wing centerline 
» ■ freestream value 


Notation 


Superscripts 
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b 
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Re 

Re 
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aerodynamic influence coefficient matrix 
aerodynamic influence coefficient 
wing span, a 

section drag coefficient > section drag/ 

wing induced drag coefficient > wing drag/ 

local skin friction coefficient , 

wing lift coefficient « wing lift/(-j P^U^C^b) 

section lift coefficient on wing centerline > 

section lift/(y P,U*cos*6C^) 

reference chord, m . 

pressure coefficient ■ (p - p„)/(j p^U^cos^B) 

minimum distance between adjacent elements 
shape factor ■ 6*/6 

co^>onent of freestream velocity normal to 
panels of kth coagionent 
static pressute, kNnm'^ 
source strength 

viscous/potential flow iteration relaxation 
factor 

Gauss-Seldel relaxation factor 
Reynolds number • U^C^j/v 

Reynolds nuaber based on momentum thickness ■ 
Ue6/v 

boundary- layer coordinate along surface in 
chordwise direction, m 
chordwise velocity, ms~^ 
spanwise velocity, ms~^ 
chordwise and spanwise wing coordinates, m 
boundary- layer coordinate normal to airfoil 
surface, m 
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^President OlkfiB ® 


I > index of viscous/potential flow iteration 
V • index of Jacobi and Gauss-Seldel iteration 

I . Introduction 


The multi-coi^onent wing is an essential element 
of the hi(h-lift system of nearly all general avia- 
tion, commercial, and military aircraft. Due to the 
high-lift requirement, the influence of viscosity 
becomes a dominant factor in determining aerodynamic 
performance. As a result, inviscid theory is not 
sufficient and viscous effects must be considered if 
accurate predictions are to be made. 

Design of these systems has generally been 
dependent on experimental verification of predicted 
aerodynamic performance. A prime example is con- 
figuration optimization where optimum slat and flap 
positions and deflections for best aerodynamic per- 
formance are generally determined thiough expensive 
and time consuming wind tunnel testing. Improved 
understanding of the aerodynamic phenomena and 
improved predictive capability would significantly 
reduce this wind tunnel effort and permit optimiza- 
tion in the early design stages. 

Stevens, Goradia, and Branden* present a theo- 
retical analysis for two-dimensional flow about 
multi-conq)onent airfoils where for the first time 
the merging of a wing wake and the upper- surface 
boundary layer of a downstream element is treated 
analytically for configurations with up to four 
elements. Bhateley and McNhirter^ and Callaghan and 
Beatty* also present analyses for two-dimensional 
multi -component configurations although neither of 
these papers treats wake-boundary layer merging 
effects. The work of Jacob and Steinbach** considers 
multi-el "'oent airfoils with flow separation and 
Ormsber jnd Chen® present a method for designing 
multi-component airfoils for maximum lift. 


1 


I 


The objective of the present work is to hej;in the 
extension of multi- component airfoil analysis into 
thrce'dimensional flows by considcrinf’ the subsonic 
flow about swept wings of infinite aspect ratio. 

The basic analytic method is outlined in Section II. 
The present analysis extends the theoretical method 
of Dvorak and Woodward^ and compares predicted 
results with experiment. The theoretical extension 
of the work of Ref. 6 includes: 1) introduction of 

advanced potential-flow solution techniques; 

2) modification of the viscous/potential flow coup- 
ling when only slight separation is present; and, 

3) use of under-relaxation of surface source distri- 
butions in the iterative coupling of the viscous and 
potential flow solutions. In the supporting experi- 
mental program, tests of a two-element configuration 
were conducted in the Ames 40- by 80-Foot Wind Tun- 
nel. The basic objective of this testin^ was to 
obtain experimental data on two-dimensional and 
infinite aspect ratio swept-wing configurations. 

The experimental program is described in Section III 
and the results are compared with the theoretical 
predictions in Section IV. 

II. Theoretical Method . 


General 


The basic theoretical method is to couple iter- 
atively potential-flow and boundary- layer analysis. 
The analysis is limited to flows with negligible 
upper-surface separation on any element. Provided 
there is no strong viscous interaction, limited 
lower surface separation closed bubbles, such as 
often occur in the cove on the lower surface of a 
wing with a slotted flap, are permitted. In the 
case of the infinite-span swept wing it is assumed 
that all spanwise gradients in the potential and 
viscous flows are equal to zero. Thus, in planes 
normal to the leading edge, the potential flow 
equations reduce to the two-dimensional form. For 
the infinite-span swept wing the chordwise boundary- 
layer characteristics are of primary interest. For 
the laminar case Jones^ shows that the cliordwise 
boundary layer is essentially independent of the 
spanwise boundary layer. For the turbulent case 
this uncoupling is not possible and the spanwise and 
chordwise boundary- layer equations must be solved 
simultaneously. 

Geometry Definition 

Dp to four elements can be analyzed, with each 
element represented by as many as 60 pairs of sur- 
face coordinates. Each individual slat or flap seg- 
ment location is related to the main wing or refer- 
ence coordinate system by pivot point coordinates 
prescribed in both the main element coordinate sys- 
tem and the individual slat or flap segment coordi- 
nate system. The selection of a rotation angle 
measured relative to the main component completes 
the specification of the element position. It is 
generally convenient to pivot a leading-edge device 
about its trailing edge and to pivot a flap about 
its leading edge although the hinge point of a flap 
could, for example, provide a ready reference point. 
Also, as part of the geometry definition, the flap 
upper-surface longitudinal radius of curvature is 
determined using cubic splines for later use in the 
finite-difference boundary-layer calculations. 


I’ot cntial Flow Analys i s 

The potential flow analysis is performed in a 
plane normal to the wing leading edge. Ihc analysis 
is a method of singularities where each clement is 
represented by a closed polygon of planer panels 
connecting the input coordinate pairs. A linear 
vorticity distribution is located along each panel 
with the requirement that the vorticity distribu- 
tion be continuous across the panel corner points. 
Thus, if there arc n panels there are n* 1 unknowns 
to be determined. The boundary condition of no flow 
through the surface, applied at each of the panel 
centers, provides n equations. The additional equa- 
tion required to close the system is supplied by 
specifying that the upper- and lower-surface velo- 
cities have a common limit at the trailing edge 
(i.e., the Kutta condition). This implies that the 
upper- and lower-surface vorticities be equal and 
opposite. For a unit freestream velocity the 
resulting system of equations to be solved for an 
airfoil with n panels is then 
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or in vector form 

(A][y1 - In] (2) 

where the nth column represents the combined influ- 
ence coefficients for the upper- and lower-surface 
trailing edge vorticity panels. 

If the trailing edge closes the vorticity must go 
to zero because the trailing edge then becomes a 
stagnation point. Although this solution is auto- 
matically given by Eq. (1) the geometry is often 
such that the influence coefficients in the nth 
column are quite small with the result that the 
matrix is poorly conditioned. In this case an 
alternate Kutta condition is used which specifies 
that the upper and lower trailing-edge vorticity 
strengths are equal to zero at the trailing edge. 

An additional unknown is supplied by introducing a 
constant strength source distribution on the surface 
of the airfoil. The resulting system of equations 
is the same as Eq. (1) except for the nth column of 
vorticity influence coefficients. That column is 
replaced by the constant source-distribution influ- 
ence coefficients and is the unknown source 
distribution strength. 


For a multi -element configuration with j compo- 
nents Eq. (2) can be written 
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self-influence coefficient tutrix 
of component i 

i )< j influence coefficient matrix of 
jth component on the ith com- 
ponent 

rjj ■ vorticity strength vector for kth component 



liquation (3) is solved by direct triungiilar 
decomposition (see, for example, Isaacson and 
Keller*) for single component airfoils. At the 
users option, triangular decomposition can also be 
used on multi -component configurations with a total 
of less than 100 panels. For multi -component con- 
figurations two iterative methods were investigated: 
block-Jacobi iteration, and block-Causs-Seidel 
iteration with relaxation. The solution is found 
by starting with some assumed initial solution 
vector. 


Y. 


(4) 


The vth iteration is then expressed in one of the 
fol lowing ways : 

1. Block Jacobi iteration: 


-(V) 

^k 


(5) 


2. Block-Gauss-Seidel iteration with relaxation: 
^k'*^ “ ♦ (1 - R) k « 1,2 j 

* * V - 1,2,... 


( 6 ) 


Equations (5) and (6) are solved using triangular 
decoaposition on the first iteration; the decom- 
posed matrices are saved for use on subsequent 
iterations. Ihe iteration process is assumed to be 
converged when 


MaxlY)* - Yj***! < e 
i«l-^ 

where e is generally taken to be 0.01. 

For configurations with two or three lifting 
elements, the Gauss-Seidel without relaxation 
(Rwl) converges in approximately one-half the num- 
ber of iterations required for the Jacobi method. 
Since the nuiBier of arithmetic operations per 
iteration are nearly equal for the two methods, the 
Gauss-Seidel method was selected as the best itera- 
tive technique. 

ORIGINAL 
OF POOR QUALOT 


As one might expect the optimum relaxation *'■ r 
for the Gauss-Scidcl method is dependent on hr 
closely the various lifting elements are coupled. 

For example, consider an airfoil with a singiv 
slotted flap. If the flap is sieved off to infinit 
the two elements are essentially uncoupled. 'Ihc 
off-diagonal block matrices of Eq. (3), wh'ch model 
the component interactions, are zero and two itera 
tions of Eq. (6) give convergence (first iteration 
gets the solution and the second is required only 
check the first). As |R-1| increases, the number 
of iterations required also increases. 

If the element; are located in close proximity, 
as is typically the case, the interaction between 
components become important. Again, consider the 
single slotted flap. Assuming the initial solution 

, k • 1,2) is a null vector, the final coupled 
solution is obtained from Eq. (6) for R«las follows’ 
1) is the solution vector for the main com- 

ponent in free air (because the upwash field fru.n 
the flap is not present, the lift level for this 
component will be less than that for the coupled 

system); 2) solution for the flap in the 

downwash field of a wing that is carrying reduced 
lift, relative to the coupled system (this reduced 
downwash field causes the flap to carry excess lift) ; 
and 3) on the next iteration the wind lift will be 
too high because of the excess upwash from the flap. 
Similarly the flap lift on the second iteration will 
be below that for the correctly coupled system. 

This oscillatory behavior indicates that under- 
relaxation should accelerate the rate of convergence 
and the amount of under- relaxation required should 
increase as the elements are moved closer together. 

An "optimum" relaxation factor given by 
R = 1 ■ exp -10 Min(I\) was determined by numerical 
i^l-Hc-l 

experiments with two and three elemeiit configura- 
tions. IS the minimum distance in fraction of 
chord between the trailing edge of the ith component 
and the upper surface of the following component. 

With this relaxation factor the number of iterntions 
required for convergence was generally reduced by a 
factor of 3 below that required with R» 1. 

Computational times required for the direct 
method of solution have been compared to the compu- 
tational times required for block-Gauss-Siedel 
iteration with relaxation method. The test cases 
were closely coupled two-component configurations. 

Ihe total number of panels used for these compari- 
sons was approximately 100. The iterative method 
obtained the converged solution in approximately 
one-half the time required by the direct method. 
Increasing the number of elements, and thus the 
total number of equations to be solved, had no 
significant impact on the number of iterations 
required by the iterative method. Thus the advan- 
tage of the Iterative metiiod over the direct will 
increase as the number of elements is increased. 

The pressure coefficient at the center of each 
panel is calculated from the surface velocity at 
that point. Tie lift and pitching-moment coeffi- 
cients are obtaineJ by integrating the pressures 
around the airfoil. 

A comparison of the present method with the 
exact conforma, mapping solution of Williams* for 


3 


■ two eleaant case is shown in Fig. 1. The angle of 
incidence is 0* and the flap deflection is 30*; 
agreesKnt with the exact solution is excellent. 



negligible influence on the overall calculation, at 
least for noderste sweep angles. Once the boundary- 
layer characteristic along the potential flow stream- 
lines have been deterained, the spanwisc and chord- 
wise paraiwters are determined. 

These chordwise boundary -layer characteristics 
are used with the correlations of Smith to deter- 
mine the point of laminar instability. These cor- 
relations relate R. to a pressure gradient parameter. 
Knowing the point of instability, Rg, and the pres- 
sure gradient; the location at which full transi- 
tion has occurred is determined by using the corre- 
lations of Granville. The Integral method of 
Cumsty and Head*** for boundary layers on infinite- 
span swept wings is used for the turbulent calcu- 
lation. 


Fig. 1 Cos;>arison of numerical and exact 
potential flow solutions: a-0*, 6^«30*. 

Boundary- Layer Analysis 

As outlined in Pig. 2 , a coii)ination of integral 
and finite difference techniques are used. Integral 
methods are used for conventional boundary layers 
because of their coqiutatlonal efficiency whereas a 
finite-difference method is used for the more com- 
plex confluent boundary- layer analysis. The stagna- 
tion line, conventional boundary- layer, and conflu- 
ent boundary- layer methods are discussed in the 
following sections. 


•ormnuLnowcMsuLaTio memeno u i n ii*— 



Fig. 2 Flow about swept, infinite-span, 
multi-element wing. 

Stagnation line Flow 

Cuq>sty and Head*** experimentally Investigated 
the stagnation- line flow characteristics of infinite- 
span swept wings. They found that the stagnation 
line boiaidary- layer integral parameters (H, 6, and 
C.) and the state (laminar or turbulent} correlate 
* j du- 

with the parameter C* ■ The stagnation 

line boundary- layer characteristics on each element 
are determined using these correlations. 

Conventional Boiaidary-Layer Methods 

Integral boundary- layer methods are used for all 
conventional bomdary layers, such as the upper- 
ond lower-surface boundary layers of all elements, 
and the upper-surface boiaidary layers of following 
elements up to the slot exists. 

The two-dimensional equations of Curie* ‘ are 
solved along external streamlines to determine the 
laminar boundary- layer characteristic. It is 
assiaed that laminar ci*oss-flow effects have a 


Should laminar separation occur prior to trans- 
ition, the correlations of Caster*^ are used to 
determine whether turbulent reattachment occurs. 

If reattachment is predicted, the calculation con- 
tinues as a turbulent flow; if catastrophic separa- 
tion is predicted, the calculation is terminated. 


Confluent Boundary- Layer Method 

The finite-difference method used to solve the 
infinite-span swept-wing, confluent boundary- layer 
equations is described by Crank and Nicholson*^ and 
by Dvorak and Head.*^ The eddy viscosity model 
used here is a modification of the two-dimensional 
method for wall Jets and turbulent boundary- layers 
developed by Dvorak.*^ These calculations include 
the effects of longitudinal surface curvature. The 
static pressure field, p(s,z), which is required 
for the solution of the boundary -layer equations, is 
determined directly from the potential flow solution. 

The initial conditions required to start the 
finite-difference calculation at the slot exit are 
constructed from: 1) the integral boundary- layer 

solution at the slot exit on the upper surface of 
the component in question; 2} the laminar potential 
core as determined from the potential-flow solution; 
and 3) the upper- and lower-surface boundary -layer 
solutions at the trailing edge of the upstream 
element. If cove separation is present the boundary 
layer at the slot exit from the lower surface of the 
upstream element is assumed to follow the one-seventh 
power law with a thickness equal to one-third the 
slot-exit width. This assumed profile is repre- 
sentative of what is observed experimentally. With 
these initial conditions and with the static pres- 
sure field known, the boundary- layer equations are 
solved in a forward marching fashion to the trailing 
edge of the component .- 


The profile drag for a streaia«ise section is 
determined by use of the method of Squire and 
Young : 


C. - 2 


H 


te 




te\“-/te 


(7) 


III. Viscous/lnviscid Coupling 

The effect of boundary- layer displacement and 
mass entrainment on the potential flow is simulated 
by distributed sources on the panels used to des- 
cribe the airfoil contour. The strengths of these 
source panels as determined direcUy from the 
boundary- layer solutions are R ■ ^ (Ug6*} where 
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u« is tha chordMlie potential-flow velocity at the 
edge of the boundary layer and i* la the boundary- 
layer dleplaceaent thlcknett given by 



A Bodlfled Kutta condition, «dilch requires that the 
flow be tangential to the trailing edge p'.iel at 
the trailing edge on the upper and lower surfaces, 
together with the boundary condition that the velo- 
city nonul to the surface be equal to the known 
source distribution, deteralnes the potential flow. 
The resulting potential-flow pressure distribution 
Is then used In subsequent boundary- layer calcula- 
tions. 

It has been observed (e.g., Brune, Rubbert, and 
Nark^*)) that conventional cyclic boundary- layer/ 
Invlsclo flow Batching usually results In a diver- 
gent or, at best, a slowly convergent Iteration 
process, especially If the viscous Interaction Is 
relatively strong (as Is often the case for hlgh- 
11ft configurations). In the present method the 
boundary- layer source distribution Is under-relaxed 
according to the formula 

q^*‘ - q* ♦ r(q* ♦ q‘) (9) 

where q^ Is the source distribution strength used in 
the 1th potential flow solution and ^ is equal to 

^ (u^S*) obtained from the ith boundary- layer 

solution. A relaxation factor of 0.3 has been found 
to be sufficient for a large variety of configura- 
tions . 


For the unswept case, separation Is defined as 
the location where the skin friction goes to zero. 
For the swept case, separation is defined as the 
point where the wall ekln-frlctlon vector forms a 
90* angle with the local potential -flow velocity 
vector. The manner in which the viscous interaction 
is handled depends on what type of separation is 
being considered. If separation occurs at some 
point, , on the upper surface of a component, 
V^°/sep 

the boundary layer calculation is necessarily ter- 
minated at that point. Based on the experimental 
observation that the static pressure is nearly con- 
stant in the separated flow, u^ is assumed constant 
and equal to Its value at the separation point. 


The displacement thickness in the separated zone 
Is obtained by linear extrapolation from the point 
of separation. Thus, in the separated zone the 


blown boundary condltlori becomes q ■ 



The method is not Intended to model extensive flow 
separation; rather, the intent Is to permit the 
calculation to continue If separation occurs, such 
as Is often encountered after the first potential 
flow solution. Separation on the lower surface of 
the last element, which seldom occurs. Is handled 
similarly. Lower-surface flow separation, such as 
might be encountered in a flap cove, is handled 
differently. Experimental observation indicates 
that even if cove separation occurs the boundary 
layer at the slot exit Is generally thin and Its 
displacement effect on the potential flow is 
expected to be minimal. Therefore, the source 
distribution in the cove-separation zone is assumed 


PArr -cj 


to be a linear interpolation made between its value 
at separation and zero at the trailing edge. 

IV. Wind Tunnel Test Program 

Model 

The rectangular planform finite wing, used in the 
experimental part of this study, is shown as it was 
mounted in the Ames 40- by 80-Foot Wind Tunnel in 
Fig. 3. This wing is equipped with a full-span, 40V 
chord, single slotted flap. The wing and flap were 
basically steel frameworks covered with wood; the 
wood in turn was coated with a glass reinforced 
plastic skin to give the desired contour. The wing 
span was 16 m with a reference chord of 1.7 m and a 
nominal extended chord of 2.1 m. The relatively 
high aspect ratio of 7.6 (based on the extended 
chord), combined with the uniform section is designed 
to provide a nearly constant span loading over the 
wing center section. In addition, the finite wing 
avoids the interference effects, caus'^d by wind- 
tunnel Side-wall junctions, that are usually encoun- 
tered with models that span the test section. These 
interference effects become especially troublesome 
at high lift and when separation plays a significant 
role. 



Fig. 3 Model mounted in 40- by 80-Foot Wind Tunnel 

The flap brackets, which remain parallel to the 
chordwise direction when the wing is yawed, permit 
horizontal and vertical flap movement. The range 
of movement in the chordwise and normal -to- chord 
direction was SV and 4V of Cq respectively. 

The basic airfoil section from which the wing- 
flap combinations were derived is an RAF J815 The 
wing-flap combinations tested arc shown in Fig. 4. 

A complete table of coordinates is listed by 
Foster et al.^* Tlic 10° flap configuration had a 
faired cove whereas the .30° flap configuration had 
no cove fairing; the leading edge of the 30° flap 
configuration was drooped 10° with the pivot point 
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on the lower surface at an of ().ir>. The 

test Kcynolds number, R^, was 3.8«10‘-'. 
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fig. 4 Airfoil configurations tested. 

Wing Forces and Moments 

Overall wing forces and moments were obtained 
fron the wind tunnel scale system and corrected for 
tunnel wall-interference effects. The total lift 
was used in conjunction with the vortex lattice 
theory of Hough^^ to determine a section angle of 
attack at the center section of the wing according 
to the formula 

a - a - L(x. ( 10 ) 

where is the induced angle of attack at the 
center section due to f ite wing effects. The 
induced angle of attack, or mean downwash at the 
wing center line, as given by the vortex lattice 
theory is 

where It Is assumed that the induced angle of attack 
scales with the measured total lift of the wing. 

Pressure Measurement : 

The center section of the wing was instrumented 
with three chordwise rows of static pressure ori- 
fices. One row was located on the model centerline, 
a second row at y ■ 0.44 Cg, and a third row at 
y ■ -0.22 Cg. No flap brackets were located between 
y = -0.22 Cg and y ■ 0.44 Cg. At each of these 
three spanwise stations there were 64 orifices on 
the main element and 30 orifices on the flap. The 
method of measuring pressures wa.s with scan! valves 
and transducers with ranges of *17kN m"^ (±2.51b/in^) ; 
the scani valves were automatically sequenced. The 
transducer output was digitized and punched onto 
data cards for subsequent data reduction on a digi- 
tal computer. 

Preliminary tests at sweep angles of 0“ and 25® 
showed no significant' spanwise gradients in the 
center section of the wing. All data presented in 
this paper were obtained on the model centerline. 

Hot-sfire Boundary- Layer Surveys 

Velocity profiles near the main element trailing 
edge and above the flap upper surface were obtained 
using a translating hot-wire probe. The drive motor 
was located inside the flap, thus minimizing the 
aerodynamic interference due to the survey device. 
Anemometer output was linearized to give a linear 
relationship between voltage and velocity and then 
recorded on a X-Y plotter as a function of probe 
position. 

Hot-wire calibrations were made using a free jet 
powered by a variable speed blower. The anemometer 


bridge wa.s temperature compensated to minimize the 
influence of wind-tunnel static-temperature varia- 
tion on anemometer output. 

V. Comparison of Theoretical and 
Experimental Resu l ts 

Section Lift 

Figures 5-7 show comparisons of the measured 
section lift characteristics, as functions of sec- 
tion angle of attack, with those determined theo- 
retically. Figure 5 is for the 10* flap configu- 
ration without sweep; Fig. 6 is for the 50“ flap 
configuration with 0* sweep; and Fig. 7 is for the 
30* flap configuration with 25* sweep. The theo- 
retical calculations have been terminated at the 
angle of attack where upper-surface separation 
first occurred. This type of separation greatly 
inhibited or prevented convergence of the viscous- 
potential flow iteration process, and resulted in 
generally unreliable results. The fact that in all 
cases separation occurred at lift levels below, but 
within 10% of, C Indicates that the existence of 
max 

separation on either the main component or flap 
provides a reasonably accurate and conservative 
estimate of maximum lift. 



Fig. 5 Section lift characteristics for 
plain leading edge: 8*0*, 6^ = 10*. 



Fig. 6 Section lift characteristics for 
drooped leading edge: 6 • 0®, 6^ * 30*. 

The general character of the results for the 
three configurations shown in Figs. 5-7 are similar. 
Comparison of the experimental measurements and the 
potential-flow solutions show strong viscous inter- 
actions resulting in observed lift coefficients well 
below those predicted by pure potential -flow theory, 
The addition of viscous interaction effects brings 
the theoretical predictions into much better agree- 
ment with the experimental results with the largest 
difference occurring, somewhat surprisingly, at 
lower angles of attack. The reason why the lift is 
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Flj{. 7 Section lift characteristics for 
drooped leading edge: 6 * 25*, = 30*. 

over-predicted at the lower angles of attack Is 
shown by the section pressure distributions pre- 
sented In Figs. 8 and 9. 


present) is suppressed; th'} result is a much 
improved agreement between theory and experiment. 
Increasing the sweep angle to 25* did not r.ienifi 
cantly alter this behavior for either the 10^ or 
30* flap configurations. 

Figure 10 compares the effect of angle of attack 
on computed integral boundary- layer characteristics 
with its effect on experimentally measured charac- 
teristics for the unswept 10* flap configuration. 
The comparisons are made for a point located at the 
trailing edge of the main component upper surface. 

The shape factor, H « ^ , is predicted accurately 

throughout the angle of attack range. Tlie tendency 
to slightly over-predlct the displacement thickness 
is caused, at least in part, bv the increased wing 
loading due to the previously mentioned cove sepa- 
ration effects . 
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Fig. 8 Section pressure distribution for 
dropped leading edge: a ■ O', 

B • 0*, 6f ■ 30*. 



Fig. 9 Section pressure distribution for 
drooped leading edge: a ■ 7.5*, 

B » 0*, • 30*. 

The flap pressure distribution for the unswept 
30* flap configuration of Fig. 8 is predicted quite 
accurately, whereas the load level for the main 
element is over-predicted. This discrepancy is 
attributed to extensive flow separation in the cove 
on the lower surface of the main component. This 
separation, althouch predicted by the boundary- 
layer analysis, results in a strong viscous inter- 
action which it not modeled accurately in the theo- 
retical analysis. 

As a result, the contribution to lift of the 
pressure distribution in the cove is over-predirted. 
In addition, the cove separation appears to have 
reduced the net circulation on the main element, 
thereby influencing the upper surface pressure 
distribution as well. At higher angles of attack, 
such as the 7.5* shown in Fig. 9, the orientation 
of the cove relative to the freestreom is such that 
the influence of cove separation (although still 
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Fig. 10 Shape factor and displacement 
thickness at wing upper surface 
trailing edge: 6 ■ 0*, 4^ • 50*. 

Velocity profiles on the flap upper surface for 
the 30* flap configuration are shown in Figs. 11 
and 12 for an angle of attack of 9*. The mean 
velocity profile in Fig. 11 is for a station 0.s‘« 
downstream of the main cosqionent trailing edge. 

The dominate feature is the wake from the main 
conq)onent. The hot-wire probe in this case did not 
get close enough to the flap upper surface to ^et 
into the very thin flap upper-surface boundary 
layer. Agreement between theory and experiment ir. 
reasonably good. Inspection of the mean velocity 
profile indicates, as assumed in the theori-t ical 
analysis, that a laminar core exits in the slot 
efflux. Also shown in Fig. 11 is the turbulence 

level , /u' ^/u^ , which provides a qualitative measure 

of the potential for turbulent transport. The 
turbulence level, significant throughout the slot 
efflux, indicates that a true laminar core is nut 
present. This turbulence is undoubtedly comprised 
of remnants of the turbulent mixing process associ- 
ated with the cove separation. The predicted and 
measured velocity profiles at the flap trailing 
edge are cosipared in Fig. 11. Although the total 
boundary- layer thickness and the point of minimum 
velocity associated with the wing wake are pre- 
dicted accurately, the experimental measurements 
show that the wing wake and flap boundary layer 
have merged to a greater extent than that predict- ' 
theoretically implying that the turbulence in thf 
so-called laminar core has a significant influencc- 
on the confluent boundary- layer development. 
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Fig. 11 Boundary- layer velocity and 
turbulence profiles downs tresn of 
wing trailing edge: a > 9*. 

B ■ 0*. • 10*. Xf/C * 0.036. 
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Fig. 12 Boundary- layer velocity profiles 
at flap trailing edge: a ■ 9*. B ■ 0*, 

(f ■ 10*. Xf/C^ * 0.39. 

Figure 13 thowi a coaparison of the conputed 
drag for the 30* flap configuration with the drag 
■aasureaents obtained from Ref. 21 for the two- 
diaansional case. Although the shape of the lift- 
drag curve is correct, the theoretical nethod con- 
sistently over-predicts the drag; this suggests 
that the Mthod of Squire and Young‘s is not 
directly applicable to multi -component airfoils. 



Fig. 13 Lift versus drag for drooped 
leading edge: B ■ 0*, • 30*. 

Figure 14 illustrates the application of the 
theoretical method to a configuration opt imitation 
problem. In this figure the effect of flap gap on 
lift at a fixed angle of attack is compared with 
the theoretical prediction for the 30* flap con- 
figuration. The optimum flap gap of 2\ of chord is 
predicted quite accurately. 



I 

<% CMOfipi 


Fig. 14 Rffect of flap gap on jcction 
lift for drooped leading edge; 
a • 6*. 8 « 0*. 

Concluding Remarks 

A theoretical method for analyting the viscous/ 
potential flow around two-dimensional and infinite- 
span multi -component airfoils has been described. 

The analysis is bas»d on iteratively coupled 
boundary- layer and potential -flow calculations. 

The viscous-flow analysis considers the coi.flwciic 
boundary l.oyer, where appropriate, unJ computes 
both spanwise and chordwise boundary- layer charac- 
teristics. In support of the theoretical program 
an experimental study was conducted in the Ames 
40- by 80-Foot Wind Tunnel. Comparisons of the 
present theoretical method with other exact solu- 
tions and with the experimental data resulted in 
the following conclusions: 

1. The present theoretical method accurately 
predicts section lift characteristics of multi- 
conq>onent configurations through moderate sweep 
angles where flow separation effects are negligible. 

2. For the configurations tested, the angle of 

attack for the first occurrence of upper surface 
separation defines a reasonably accurate and con- 
servative estimate of C at sweep angles of 0* 
and 25*. ^max 

3. The present method accurately predicts the 
optimum flap gap for maximum lift at a fixed angle 
of attack. 

4. In addition to the study of two-dimensional 
and infinite-span swept multi-coBq>onent airfoils, 
this method will be useful in the analysi. .>> t!., 
high-lift characteristics of relatively higl. aspi.t 
ratio, moderately swept, finite wings. 

5. The influence of cove separatioi. o,. .. 

state (lamina’* or turbulent) of the slot eiiiux can 
have a strong influence on the development of the 
confluent boundary layer over the flap upper sur- 
face. 

6. The method of Squire and Young for computi” ' 
drag does not (in its original form) accuratel; 
predict the drag of multi -component configurati 

7. Potential -flow matrix computer solution 
times for multi -component configurations are sig- 
nificantly reduced, from those required for direct 
trianeular decomposition solution, by using blocV. 
CausS'Siedel iteration with relaxation. 

8. The present potential- flow method, wVi’r 
utilizes linear vorticity, gives excellent agree- 
ment with exact solutions for the inviscid flow 
around multi-cosqwnent configurations. 
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